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Abstract 



A new next-to-leading order Monte Carlo program for the calculation of 
jet cross sections in photoproduction is described. The contributions from 
both resolved and direct components are included to 0{aa1). Properties of 
the predictions for various inclusive jet and dijet observables are discussed 
and comparisons with HERA data are presented. 

PACS number(s): 12.38.Bx,13.60.Hb,14.70.Bh 
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I. INTRODUCTION 



Electromagnetic interactions have long been used to study both hadronic structure and 
strong interaction dynamics. Examples include deep inelastic lepton-nucleon scattering, 
hadronic production of lepton pairs, the production of photons with large transverse mo- 
menta, and various photoproduction processes involving the scattering of real or very low 
mass virtual photons from hadrons. In particular, the photoproduction of jets with large 
transverse momenta is calculable in QCD and offers additional complementary information 
to that obtained from the study of the hadroproduction of jets In the photoproduction 
case there are contributions where the photon's energy contributes entirely to the hard- 
scattering subprocess; these are often referred to collectively as the direct component. In 
addition, a real photon can interact via its hadronic substructure. These contributions com- 
prise the resolved component, a review of which is contained in Ref. p[. Therefore, the 
photoproduction of jets allows one to investigate new production mechanisms, probe the 
hadronic substructure of the photon, and study the conventional hadronic jet production 
mechanisms which these processes have in common with the hadroproduction case. 

Two basic approaches are commonly employed for generating predictions for hard- 
scattering processes. On the one hand, it is possible to perform calculations for a specific 
observable in which the integrals over the subprocess variables are done analytically, leaving 
only the convolutions with the parton distributions to be done numerically. On the other 
hand, if the subprocess integrations and the parton distribution convolutions are done using 
Monte Carlo techniques, it is possible to generate predictions for a variety of different ob- 
servables simultaneously. This latter approach is sometimes referred to as a fully differential 
Monte Carlo calculation. Several groups have performed next-to-leading-order calculations 
of jet photoproduction in varying degrees of generality using one or the other of these two 
approaches. In 0, [^], and subprocesses which involved the photon were kept up to 
0{aal). Thus, the direct component was calculated to next-to-leading-logarithm accuracy 
while the resolved component was calculated in the leading-logarithm approximation. The 
resolved component was calculated to next-to-leading-logarithm accuracy in , , |^ , and 
for single inclusive production. For a review see l]!^ . Recently, both direct and resolved 



contributions calculated to next-to-leading-logarithm accuracy have begun to appear [jTT 
|T^, [0 in a fully differential form. 

The purpose of this work is to present a discussion of a calculation which is based on the 



phase space slicing method using two cutoffs [14|. Both the direct and resolved components 



are included to next-to-leading-logarithm accuracy. The result is fully differential and im- 
plemented in a Monte Carlo style program which allows the simultaneous histogramming of 
many distributions incorporating experimental cuts. It represents an improvement of earlier 
results which included the direct component only at NLO and an elaboration of the very 
brief results for both components already presented [|TT|. Details of the calculational method 



are presented as well as comparisons with recent data. Some comments on various unsettled 
issues that arise when comparing with dijet cross sections are also given. 

The remainder of the paper is as follows. The phase space slicing method is reviewed in 
Sec. II. Numerical results are compared with HI and ZEUS data and related physics issues 
are discussed in Sec. Ill, while the conclusions are given in Sec. IV. 
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II. METHOD 



In this section we describe the calculation of QCD corrections to two-jet production 
in electron-proton scattering using the phase space slicing method. Before discussing the 
QCD corrections, it is necessary to recount the connection between electron-proton and 
photon-proton scattering. For small photon virtualities the two are related using the 
Weizsacker- Williams approximation wherein one assumes that the incoming electron 
beam is equivalent to a broad-band photon beam. The cross section for electron-proton 
scattering is then given as a convolution of the photon distribution in an electron and the 
photon-proton cross section 
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The integration limits and the maximum photon virtuality Q^ax determined from the 
(anti) tagging conditions of the experiment. The energy of the photon is related to the 
energy of the incident electron by = yE^. We used the improved photon distribution in 
an electron [|l^ given by the formula 
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where rrie is the electron mass and ttgm is the electromagnetic coupling. Within this approx- 
imation, QCD corrections to electron-proton scattering correspond to QCD corrections to 
photon-proton scattering to which we now turn. 

In this version of the phase space slicing method 



IM two small cutoffs 5. and 5r are 



used to delineate regions of phase space where soft and coUinear singularities occur. Let the 
four vectors of the three-body subprocesses be labeled pi + P2 ^ + P4 + Ps, and define 
the Mandelstam invariants Sij = {pi + pjY and tij = {pi — pjY- Consider, for example, 
the ■j{pi)g{p2) Q'(P3)^(P4)5'(P5) subprocess whose matrix element becomes singular as the 
energy E^ of the final state gluon becomes soft. Define the soft region S hj < E^ < 
^sy/suf^ and the complementary hard region H by 5sa/si2/2 < E^ < y/suf^, both in 
the pi + p2 rest frame. The two-to-three body contribution to the cross section is then 
decomposed as 
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where |Mp is the three-body squared matrix element and dT^ is the three-body phase space. 
Within S one sets P5 = everywhere except in the denominators of the matrix elements 
and then analytically integrates over the unobserved degrees of freedom in n space-time 
dimensions. The result, proportional to the leading order cross section, contains double and 
single poles in n — 4, and double and single logarithms in the soft cutoff 6s- Terms of order 
6s are neglected. Next, the collinear regions of phase space are defined to be those where 
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any invariant {sij or tij) becomes smaller in magnitude than ScS^. The hard region is then 
decomposed into collinear C, and non-coUinear C, regions as follows: 



2s 12 Jh 



JH 2Si2 JHC 2Si2 JHC 



(2.4) 



Within HC one retains only the leading pole of the vanishing invariant in the squared 
matrix elements. Exact collinear kinematics are used to define the integration domain of 
HC, which is valid so long as 5c <C Sg. The integrations over the unobserved degrees 
of freedom are then performed analytically in n space-time dimensions giving a factorized 
result where single poles in n — 4, and single logarithms in both cutoffs 5c and 5^, multiply 
splitting functions and lower-order squared matrix elements. Terms of order 5c and 5s 
are neglected. The soft and final state hard collinear singularities cancel upon addition 
of the interference of the leading order diagrams with the renormalized one-loop virtual 
diagrams. The remaining initial state collinear singularities are factorized and absorbed 
into the parton distributions. The integrations over the singularity-free portion of the three- 
body phase space HC are performed using standard Monte Carlo techniques. When all of 
the contributions are combined at the histogramming stage the cutoff dependences cancel, 
provided one looks at a suitably defined infrared-safe observable. An example will be given 
below. 

The matrix elements squared for all two-to-two and two-to-three parton-parton scattering 
subprocesses through 0{al) are from the paper of Ellis and Sexton and those for the 
photon-parton scattering subprocesses through 0{aa1) are from the paper of Aurenche et 

ai. m- 

In order to compare with experimental data, a suitable jet definition must be chosen. 
For this study, the Snowmass algorithm [jl9| has been used, supplemented with an i?sep 
constraint |^ as follows. Let rji (rjj) be the pseudorapidity and 4>i be the azimuthal 
angle of parton i (j) and then define 



Ri 



Partons i and j are merged if 
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Rij < min 
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Note that i?sep = 2i? corresponds to having no -Rsep constraint. The 4- vector recombination 
is of the type 



E^^ = Ej^ + Ej^ 

r/j"* = {r]'E'j, + 7]^E^T)/E^T 
= (0*E^ + 0JE^)/Ej?\ 



(2.7) 



Other types of jet algorithms may also be used. 
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FIG. 1. The single jet inclusive cross section integrated over 30GeV < E}^^ < 40GeV and 
< rf^^ < 2. The two-body (negative) and three-body (positive) contributions together with their 
sum are shown as a function of 6s with 6c = 6s /a for a = 100 (dash line) and a = 200 (solid line). 
The bottom enlargement shows the sum for a = 100 (circles) and a = 200 (diamonds) relative to 
±5% (dot line) of the average of the last four points (solid line). 

III. RESULTS 
A. General properties 

Using the results of the method described in the previous section we have constructed 
a program to calculate jet photoproduction cross sections. The program uses Monte Carlo 
integration so it is possible to implement experimental cuts, provided that they are defined 
in terms of partonic variables. 

The CTEQ4 E^ proton parton distribution set was used throughout, in conjunction 



with either the GRV p2| or GS photon set. Unless otherwise noted, the GRV set 
was used. The calculation was performed in the MS factorization scheme. Therefore, the 
corresponding MS distribution sets were used when working at next-to-leading order (NLO). 
For leading order (LO) results, the appropriate LO sets were used. Similarly, the two-loop 
version of the strong coupling as was used with matching across quark thresholds for the 
NLO results, and at LO the one-loop value was used. The value of A^^° was taken from 
the proton parton distribution set. 

For future use we define three sets of beam energy and tagging conditions as follows: 

1993 HI where = 26.7 GeV, Ep = 820 GeV, Q^^^ = 0.01 GeV^ and 0.25 <y < 0.7, 1993 
ZEUS where E^ = 26.7 GeV, Ep = 820 GeV, Q^^^ = 4 GeV^ and 0.2 <y < 0.85, and finally 

1994 ZEUS where E^ = 27.5 GeV, Ep = 820 GeV, Q^x = 4 GeV^ and 0.25 <y< 0.8. 

As explained above, when contributions from the two-body and three-body pieces are 
combined to form a suitably defined infrared-safe observable the cutoff dependences cancel. 
This provides a check on the method and also on its implementation. As stated in Sec. 
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Ef^(GeV) 



FIG. 2. Scale dependence of the single jet inclusive cross section as a function of integrated 
over — 1 < < 2. The LO (dash line) and NLO (solid line) results are shown for two different 
scale choices, /x = (top) and fj, = 2E^^ (bottom). 



II, the method requires that dc ^ ^s- Therefore, we take Sc = Sg/a with a ~ (9(100) and 
verify that the result is independent of a. As an example, the single jet inclusive cross 
section integrated over 30 GeV < Ej^^ < 40 GeV and < 7]^"^^ < 2 is shown in FIG. 1 for the 
1994 ZEUS conditions with R — 1 and i?sep = 2. The two-body (negative) and three-body 
(positive) contributions together with their sum are shown as a function of 5s for a = 100 
(dash line) and a = 200 (solid line). The bottom enlargement shows the sum for a = 100 
(circles) and a = 200 (diamonds) relative to ±5% (dot line) of the actual result (solid line). 
For this particular observable, the result is sufficiently independent of the cutoffs below 
5s = 10-2. 

One of the primary reasons for carrying out a NLO calculation is that of reducing fac- 
torization and renormalization scale dependence. This is illustrated in FIG. 2 for the single 
jet inclusive cross section as a function of integrated over — 1 < r]^^^ < 2 for the 1993 
HI conditions with R = 1 and i?sep = 2. The factorization and renormalization scales have 
been set equal to The LO (dash line) and NLO (sohd line) results are shown for two 
different scale choices // = (top) and = 2E^^ (bottom) where E^^ is the maxi- 

mum transverse energy associated with the rel event two- or three-body weight. The single 
jet inclusive cross section as a function of rj^^^ integrated over i?^* > 15 GeV is shown in 
FIG. 3. In both cases, the NLO calculation shows reduced scale dependence relative to the 
LO result. 

As stated in the introduction, it is conventional to describe jet photoproduction in terms 
of two components. The direct or point-like component corresponds to the case where 
the photon participates wholly in the hard-scattering subprocess, whereas the resolved or 
hadron-like component corresponds to a situation where the photon interacts as if it con- 
tained partons. In the direct component there is, then, no remnant beam jet and this 
leads to a simpler event structure than in the purely hadronic case. The decomposition is. 
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FIG. 3. Scale dependence of the single jet inclusive cross section as a function of rf'^'^ integrated 
over ^Jj?* > 15GeV. The LO (dash line) and NLO (solid line) results are shown for two different 
scale choices n = ET^^^ (top) and \jl = 2Ejl^'^ (bottom). 




FIG. 4. The NLO single jet inclusive cross section as a function of integrated over 
— 1 < < 2 in three different regions < x-y < 1 (solid), > 0.75 (dot), and Xj < 0.75 
(dash). 
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FIG. 5. The NLO single jet inclusive cross section as a function of tt^''* integrated over > 17 
in three different regions < < 1 (solid), > 0.75 (dot), and Xj < 0.75 (dash). 




Mjj (GeV) 

FIG. 6. Dijet cross section as a function of the invariant mass of the two highest Et jets 
decomposed into the regions < x-y < 1 (solid), x^ > 0.75 (dot), and x^ < 0.75 (dash). 
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however, an artifact of viewing the process from the perspective of leading order QCD. At 
next-to-leading order , the two components mix and a complete separation is not possible. 

Although a strict separation into direct and resolved components is not possible, one 
can still define kinematic regions where either component dominates. One such definition 



2M uses the observed momentum fraction of the parton coming from the photon which 



is given by 

= (^j^^ie-^J^'i + Ej.'*2e-''J^*2)/2E^. (3.1) 



Here Et and t] are those of the two highest transverse energy jets. The region > 0.75 
corresponds primarily to the direct contribution and < 0.75 to the resolved contribution. 

It is interesting to see how the two components populate phase space. To this end, 
shown in FIG. 4 is the NLO single jet inclusive cross section as a function of Ej^^ integrated 
over —1 < if^^ < 2 in three different regions < < 1 (solid), x^ > 0.75 (dot), and 
x-y < 0.75 (dash) for the 1993 ZEUS conditions with R = 1, and i?sep = 2. The cross 
section dominated by the resolved component has a steeper Et dependence than that which 
is dominated by the direct component. This is due to the extra convolution associated 
with the photon parton distribution. In other words, only a portion of the photon's energy 
contributes to the hard scattering for the resolved component. Hence, at large Et values 
the direct component dominates. In FIG. 5 the NLO single jet inclusive cross section is 
shown as a function of rj^^^ integrated over Et^ > 17GeV. The direct component dominates 
at negative rapidity whereas the resolved component takes over in the positive direction. 
In the coordinate system used in the HERA experiments, positive rapidity corresponds to 
the direction of the proton beam. Therefore, if the observed jet is in the positive rapidity 
region, this corresponds to larger x values from the proton and smaller ones from the photon. 
Conversely, if the jet is in the negative rapidity region, larger photon x values and smaller 
proton X values are favored. Since the contribution dominated by the direct component has 
x^ > 0.75, it dominates in the negative rapidity region, as shown. In FIG. 6 the dijet cross 
section is shown as a function of the invariant mass of the two highest Et jets. Again, the 
steeper resolved component dominates over the flatter direct component at low Mjj. 



B. Comparison to single inclusive data 



Having discussed the general properties of the predictions, it is of interest to see how 
they compare with existing data. The first comparison is to the single jet inclusive cross 
sections as measured by HI |^. The single jet inclusive cross section as a function of i?^* 
along with our NLO result (GS photon: solid lines, GRV photon: dash lines) for the 1993 
HI conditions with fi = Et^^, R = 1, and Rsep = 2 is shown in FIG. 7. The data and theory 
have different slopes with the theory falling below the data at low Et^ for both photon sets. 
The agreement is slightly better for the case where t]-'^* < 1. This is seen more clearly in 
FIG. 8 which shows the single jet inclusive cross section as a function of V''* for ^j?* > E:^^ 
with Et^^ = 7, 11, and 15 GeV. In the proton direction the theory consistently undershoots 
the data. The agreement improves as i?™™ increases. One possible explanation for the 
observed discrepency may be the need to increase the gluonic content of the photon which is 
poorly constrained in this x region. However, this would be premature as underlying event 
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FIG. 7. The single jet inclusive cross section as a function of E'j'* ^ measured by HI |25 
compared with our NLO result (GS photon: solid lines, GRV photon: dash lines). 
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FIG. 8. The single jet inclusive cross section as a function of if'^^ as measured by HI |25] 
compared with our NLO result (GS photon: solid lines, GRV photon: dash lines). 
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contributions, which were not removed from the data set or included in our calculation, may 
also provide an explaination. 

We next switch to 1993 ZEUS conditions and compare to the published ZEUS data 
on single jet inclusive cross sections. In FIG. 9 the measured single jet inclusive cross 



section is shown as a function of £"5?* along with our NLO result (GS photon: solid lines, 
GRV photon: dash lines) with fi = E^^^, R = 1, and Rsep = 2. In this case the agreement 
with the data is somewhat better than that shown in FIG. 7. This is further illustrated by 
the comparison with the single jet inclusive cross section as a function of rf'^^, shown in FIG. 
10. The trend, although less pronounced, is the same as with the HI data. 



C. Comparison to dijet data 



Data for several different types of observables involving two large transverse momentum 
jets have been published. In the ZEUS dijet analysis a sample of dijet events is obtained 
by requiring i?;^*^ > i?™". In principle, this presents a problem in the following 

sense. The pair of jets can be accompanied by additional soft gluons. The singularities 
associated with the soft gluon emission cancel infrared singularities coming from one-loop 
contributions. This cancellation has already taken place in the present calculation. That is, 
the soft gluons with energies up to 6sy/si2/2 have already been integrated out. After the 
cancellation there are residual In 6s terms in the two-body portion of the calculation which 
cancel against logarithms built up by integrating over the unobserved gluon radiation in 
the three body contribution. The requirement that both jets have Et > restricts a 

portion of this integration and leads to an incomplete cancellation of the In 6s dependence. 
We have, nevertheless, used this definition for the dijet sample, having found that the cutoff 
dependence of the results is negligible. However, the problem still exists, at least in principle. 
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FIG. 10. The single jet inclusive cross section as a function of rf^^ as measured by ZEUS |26] 
compared with our NLO result (GS photon: solid lines, GRV photon: dash lines). 
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FIG. 11. Dijet cross section for > 0.75 as a function of ry integrated over 
for i?™™ = 6,8, 11, 15 GeV as measured by ZEUS The solid curves correspond to requiring 

> ^min ^^^i^ ^i^g dashed curves correspond to the cuts used in |l^, as described in 

text. 
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FIG. 12. Dijet cross section for 0.3 < a^-y < 0.75 as a function of r] integrated over Ej^^ > E'™™ 
6, 8, ll,15GeV as measured by ZEUS [p7|]. The solid curves correspond to requiring 



for = _ 

E^^^\Ei^^^ > ^^in ^ijiig ^i^g dashed curves correspond to the cuts used in |T3], as described in 
text. 
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< 0.5 are shown in FIGS. 11 and 12 for e}^^\E^^^ 



> Ef'' with 

6, 8, 11, and 15 GeV. The solid curves are our results for the 1994 ZEUS conditions 



using the GRV photon set, /i 



TTTnax 



and R = 1. A value of R 



1 was used following 

28| where it was found that this value used at the parton level most closely corresponds 



Sep 



to the KTCLUS algorithm used to define the jets. No systematic cutoff dependence was 
observed for these results, within the limitations of the Monte Carlo errors. The solid curves 
display the same trend that was observed in the single inclusive plots discussed previously 
insofar as the direct component comes closer to the data than does the resolved, the data are 
underestimated in the low Et region, and the agreement gets better as ii^™™ increases. This 
is consistent with there being some residual underlying event contribution to the observed 
jets which has not been removed during the data analysis. 

The problem of defining a sample of two jet events in a manner which is amenable to 



a next-to- leading order theoretical treatment has been studied in [|T^]. In order to avoid a 
dependence on the theoretical cutoffs used in the course of the calculation, it was proposed 
to calculate the dijet cross section by allowing the second jet to have a transverse energy less 
than E'™™ if the third unobserved jet is soft, i.e. has a transverse energy of less that 1 GeV 
(i?;^*^ < IGeV). This approach slightly expands the available phase space when both jets 
have transverse energies near ii^™™ and, therefore, avoids the problem of a potential cutoff 
dependence. Results using this algorithm are shown in FIGS. 11 and 12 by the dashed lines. 
The relaxation of the i?^*" constraint on the second jet generates additional contributions 
from the 2 — >■ 3 subprocesses so that the dashed curves lie above the solid ones and are in 
better agreement with the data. However, the cuts do not correspond to those utilized in 
the definition of the data sample. Given the unresolved question of the underlying event 
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FIG. 13. Dijet angular distribution da/d] cos9*\ as measured by ZEUS 
at COS0* = compared with LO (dash hnes) and NLO result (solid lines). 
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contributions, it is premature to draw detailed conclusions from comparisons with these 
data. 

Angular distributions in dijet production are sensitive to the spin of the exchanged parti- 
cle. Let 9* be the scattering angle of the two highest Et jets in their center of mass. In LO, 
the direct component is dominated by t— channel quark exchange which has a characteristic 
angular dependence given by (1 — | cos^*|)~^. The resolved component, however, is dom- 
inated by t— channel gluon exchange and has a (1 — | cos6'*|)~^ angular dependence. This 
property is preserved at NLO as demonstrated by the solid lines in FIG. 13 where the dijet 
angular distribution da/d\ cos^*| is shown as a function of | cos6'*| in the regions > 0.75 
(left) and < 0.75 (right) for the 1994 ZEUS conditions with R = 1, and i?sep = 2. Also 



shown in the figure is the ZEUS ||2^ measurement and the LO (dash line) result. To better 
compare shapes, the distributions are normalized to one at |cos^*| = 0. The shapes are 
very insensitive to how one handles the dijet definition issue discussed above. The excel- 
lent agreement between the predictions and the data confirms the theoretically expected 
behavior. 



IV. CONCLUSION 

In this paper, results from a new next-to-leading-order program for jet photoproduction 
have been presented. The calculation was performed using the phase space slicing method 
and implemented in a Monte Carlo style program. Physical observables are independent of 
the phase space slicing parameters. 

As is the case in other jet production processes, the renormalization and factorization 
scale dependence is reduced relative to leading-order predictions. A decomposition into 
direct and resolved components was also studied using a definition based on the observed 
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momentum fraction of the parton from the photon. The direct component dominates at 
large Et and in the photon direction while the resolved dominates at low Et and in the 
proton direction. 

Comparisons with single inclusive jet cross section measurements by HI and ZEUS as a 
function of E^^^ and rf^^ show that the theory tends to lie somewhat below the data in the 
forward rapidity region. This is the region corresponding to low values for and which 
is dominated by the resolved component. One might suspect that this indicates a need for 
increasing the size of the photon parton distributions in this x range, but this conclusion 
would be premature since contributions from underlying event structure could account for 
this, as well. 

Our results were also compared with dijet cross section measurements from ZEUS. The 
cuts used to define the dijet sample are such as to restrict a portion of the phase space 
needed to insure the proper cancellation of the cutoff dependence in the two- and three-body 
contributions. In practice, however, the variation was found to be smaller than the statistical 
errors associated with the Monte Carlo integrations. This point aside, the predictions for 
the normalized angular distributions agree well with the data. The 77 distributions show the 
same pattern of deviations from the theoretical predictions as was observed in the single 
inclusive jet measurements. 

As additional data are acquired it is possible that jet photoproduction will be used to 
place additional constraints on the parton distributions in the photon. Such data would 
complement data for the photon structure function which are primarily sensitive to the 
quark distributions in the photon. The jet photoproduction process, on the other hand, 
directly involves the gluon distribution in the photon. However, in order for such data to 
provide useful constraints on the photon parton distributions, the effects of the underlying 
event structure must be minimized. Higher Et thresholds combined with a smaller cone size 
in the jet definition would help in this regard. 
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